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A search for radio pulsars around low-mass white dwarfs 
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ABSTRACT 

Low-mass white dwarfs can either be produced in low-mass X-ray binaries by stable 
mass transfer to a neutron star, or in a common-envelope phase with a heavier white 
dwarf companion. We have searched 8 low-mass white dwarf candidates recently iden- 
tified in the Sloan Digital Sky Survey for radio pulsations from pulsar companions, 
using the Green Bank Telescope at 340MHz. We have found no pulsations down to 
flux densities of 0.6-0.8 mJy kpc -2 and conclude that a given low-mass helium-core 
white dwarf has a probability of < 0.18 ±0.05 of being in a binary with a radio pulsar. 
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1 INTRODUCTION 

High-precision timing of millisecond and binary pulsars has 
applications ranging from setting li mits on the stochasti c 
gravitational- wave background (e.g. IJaffe fc Backer! [2003). 
to studying binary evolution, neutron-star mas ses and equa- 
tions of state (e.g. ILattimer fc Pra kash 200 l]), and testing 
the predictions of gene ral relativity fe.g. lTavlor fc Weisberd 
19891: iLvne et ~al l2004h . Wide-area systematic searches have 
historica lly been very succes sful in finding millisecond pul- 
sars fe.g. lLorimer et all 1990 ). In contrast to these wide-area 
surveys, directed pulsar searches use typical associations of 
neutron stars with other objects to increase the chance of 
finding a pulsar in a given telescope pointing. This allows di- 
rect searches to use fewer pointings, with longer integration 
times and hence better sensitivity. Directed searches towards 
certain object types also elucidate their evolutionary rela- 
tion wit h radio pulsars through the i ncidence of the associa- 
tions (cf. lStaelin fc R cifcnstcin 1968]). The classic example is 
searching for young pulsars toward pulsar wind nebulae and 
supernova re mnants, often with known central X-ray point 
sources (e.g. ICamilo et al]|2002D . and the success of these 
searches has enabled meaningful statistics on the interac- 
tions between the ra dio pulsars and the supernova remnants 
fe.g. lChevaiieill2005th Directed searches for recycled millisec- 
ond pulsars have so- far been limited to globular clusters (e.g. 
iRansom et alj|2005l) and a few isolated steep-spe ctrum tar- 
gets fe.g. lBacker et al.|[l982l : iNavarro et al.lll995l ). 



The situation has changed with the publication of 
cat alogs of white dwarfs from deep surveys. In particu- 
lar, [eg cmman et al.l l|2004 ) list over 2500 spectroscopically 
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Figure 1. Measured temperatures and surface gravities plot- 
ted onto simulated evolutionary tracks for low-mass helium-core 
white dwarfs of different mass. Shown are the 7 candidates (open 
circles) and 2 pulsar companions (filled squares) that fal l withi n 
the simulated range. Data taken fro m Klcinm an et al ] l|2004l) . 
Eiscnstcin ct al. I j2006h andlBassa et al.l ||2006|) , pl otted over evo- 
lution ary tracks from lAlthaus et al. ~ll200ll) and ISerenelli et all 
l|200lh as cited in lLiebert et al] l|2004h . 



identified white dwarfs from th e first data relea se of the 
Sloan Digital Sky Survey (SDSS. lYork et al.ll2000h . Among 
these are a handful which appear to have very low surface 
gravities (log g ^ 7.5) and fairly high effective tempera- 
tures of 2-3xl0 4 K, indicating probable masses of less than 
0.45 Mq and a helium composition ( Figure [p. These ob - 
jects are described in more detail in iLiebert et al.l l|2004h . 
in which one object, SDSS J123410.37-022802.9, is shown 
to have an excellently determined spectrum and a likely 
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Table 1. Upper limits were derived using the longest individual 
integ ration time, a 22 K receiver temperature, sky temperatures 
from lHaslam et al. I jl98lT) using a spectral index of —2.6, a min- 
imum detection SNR of 10, a 2.0 K Jy — 1 telescope gain and a 
constant 10% pulsar duty cycle. 



mass of ~ 0.18-0.19 Mq. Such a low-mass, helium-core 
white dwarf is the prototypical companion of a highly re- 
cycled millisecond pulsar: in the low-mass X-ray binary 
phase much of the white-dwarf envelope is transferred to 
the neutron star, spinning it up. While this is not the 
only plausible explanation for the existence of this star (a 
white-dwarf companion is also a potential explanation), the 
match with the masses of known pulsar com panions such as 
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.van Straten et al. 
Bassa et alj|2003h and 
20061 ) is tan 



those of PSRs J1012+5307 (0.16 ± 0.02A/ Q 
19981) . J0437-4715 (0.236 ± 0.017^ 
200 lfi . J0218+4232 (0.20 ± 0.02M« 
J1911-5958A (0.175 ± O.Q1OM , bassa et al 
talizing. 

Using the Green Bank Telescope (GBT) we have 
searched for radio pulsations from possible neutron-star 
compa nions to the 8 white dwarfs listed in iLiebert et al.l 
(|2004h and Tabled We describe our observing setup in § [51 
our results in § [3] and discuss our findings in § [4] 



2 OBSERVING SETUP 

All sources were observed with the GBT on two or more 
occasions between October 21 and 26, 2004. At a cen- 
tral observing frequ ency of 340 MHz the Spigot backend 
l|Kaplan et al.l 120051 ) provided 50 MHz of bandwidth di- 
vided into 1024 spectral channels, recording total inten- 
sity in 16 bits every 81.92 /is. Sessions started with 1- 
minute observations of bright test pulsar PSR J1012+5307 
or PSR J2317+1439, followed by several 23- minute (2 24 sam- 
ples) or 46-minute (2 25 samples) searches in the direction of 
the white dwarfs listed in Table Q] 

Over the following mont hs the data we re processed us- 
ing the SIGPROC package (lLorimerl l200ll ) on the GASP 
computer cluster IjDemorest et al.l I2004T ). As described in 
more detail below, data were examined for radio interfer- 
ence, corrected for interstellar dispersion and binary accel- 
eration, and searched for periodic signals. 

Each observation was checked for interference by sum- 
ming all frequency channels at zero dispersion measure 
(DM), searching for the strongest periodicities in the result- 
ing time series and then localizing the matching interference 
within the data. Usually several affected frequency channels 
or data segments had to be removed, mostly containing non- 
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Figure 2. Detection signal-to-noise versus trial dispersion mea- 
sure and acceleration for a partly simulated 40-minute observa- 
tion of PSR J1 012+5307. U s ing th e PSR J1012+5307 binary pa- 
rameters from lLange et al,l 11200 If ), we extrapolated one of our 
1-minute test observations to a 40-minute observation around or- 
bital phase 0.25. These data were subsequently processed in our 
search pipeline which detected the pulsar at an acceleration of 
4.8ms -2 and a DM of 9.0pcc m~ 3 , close to the literature value 
of 9.0233 llNicastro et alj|l995h . 



periodic bursts or (harmonics of) 1.266 and 60 Hz interfer- 
ence. The source of the former frequency is unclear. 

We next formed multiple trial time series for DMs be- 
tween 5pccm -3 and twice the maximum expected DM in 
the directio n of each white dwarf, w hich is on average around 
40pccm" 3 l|Cordes fc Laziol 120021 ) . using 0.1 to 0.5pccm~ 3 
intervals. 



We are specifically targeting millisecond pulsars in bi- 
nary systems and the acceleration that accompanies the 
pulsar's orbital motion can introduce significant changes in 
the apparent pulse period, hampering detection of a pe- 
riodic pulse signal. The maximum acceleration of 90% of 
all currently kno wn pulsars falls within the —25 ms -2 to 
+25 ms -2 range ((Manchester et aT]l2005l ) and we therefore 
resample each dedispersed time series for 250 trial accelera- 
tions throughout this range. 

For each dedispersed and resampled time series we then 
compute the discrete Fourier transform and look for peaks, 
caused by periodic signals. As different pulse shapes and 
duty cycles produce peaks at different harmonics of the pulse 
frequency, we also fold the spectrum 2,4,8 and 16 times to 
add these harmonics. For each peak in these 5 raw and folded 
spectra that has a signal-to-noise ratio (SNR) exceeding 5 cr, 
the signal frequency and SNR are stored. 

After candidates have been identified in this fashion for 
all trial dedispersion and acceleration values, occurrences 
of the same period at different DMs and accelerations are 
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Figure 3. Estimated minimum detectable flux S m i„ as a func- 
tion of period, for 4 characteristic dispersion measures (0, 20, 40 
and 60pccm -3 ). Curves are calculated for the observing setup 
described in § [2] and Table [T] for a 2800 s integration time and 
30 K sky temperature. 



combined to summarize and plot the signal properties. The 
first diagnostic is a two dimensional plot of SNR versus 
DM and acceleration, as shown in Figure [2] for test pul- 
sar PSR J1012+5307. In order of peak SNR these plots are 
inspected for a well-defined response in both DM and accel- 
eration. For the second diagnostic the time series are dedis- 
persed, de-accelerated and folded at the best values for DM, 
acceleration and period. This allows us to check on pulse 
shapes, signal bandwidth and signal transience or persis- 
tence. 



itor of the observed white dwarf interacts with a close bi- 
nary companion and loses its envelope. In one channel of 
evolution, the progenitor's companion is a heavier white 
dwarf and a common-envelope phase exposes the progeni- 
tor's helium core. As the companion white dwarf is older 
and cooler it is much dimmer and harder to detect than 
the observed white dwarf. In the other channel of evolution 
the progenitor of the observed white dwarf (WD) is in a 
binary with a neutron star (NS). As the progenitor evolves 
onto the red giant branch it overflows its Roche lobe and a 
phase of mass transfer starts that both recycles the pulsar 
and removes part of the progenitor's h ydrogen envelope (cf. 
iBhattacharva fe van den Heuvellll99l[ ). 

We will outline the factors that determine radio pul- 
sar detectability to estimate the probability Pmsp that a 
individual helium-core white dwarf was formed in the pres- 
ence of a neutron star, as opposed to formation near another 
white dwarf. 



We first evaluate our luminosity completeness. The 
1 " 0( ' combination of our radio-flux upper limits with the white- 
dwarf distance estimates allows us to constrain the radio- 
pulsar luminosities. Typically, our candidate white dwarfs 
are 3-4 magnitudes brighter in r than t he companion to 
PSR J0218+4232 ( 18.2- 20.5 versus 23.9, lEisenstein et all 
120061 ; iBassa et al] I2003T I and hotter by a factor of 2. 
We therefore use the low e nd of the distan ce range to 
PSR J0218+4232 (2.5-4 kpc, iBassa et al]|2003h to estimate 
the distance for each candidate to be 2 kpc. At this distance 
our radio-flux upper limits amount to a luminosity of 0.6- 
0.8 mjy kpc -2 . As 95% of curre ntly known field millisec ond 
pulsars are brighter than this (jManchester et alll2005l ) we 
assume we are complete to a level of P; um =0.95. 

Millisecond pulsars generally have relatively wide 
beams that show up as high duty cycles (~20% on aver- 
age) and increase the chance that the beam sweeps over the 
Earth. For millisecond pulsars this beaming fraction Pb ea m 
is thought to be as high as 0.7±0.2 IjKramer et al.l l 1998). 



3 RESULTS 

None of the observations were found to contain signals that 
merited follow-up. In Table Q] we estimate the pulsar-flux 
upper limits for the longest observation of each candidate. 
The upper limits quoted are for pulsars with periods of 20 ms 
or more, for which our search algorithm is most sensitive; in 
Figure [3] we show the typical dependence of this minimum 
detectable flux S m i n on the pulsar period and dispersion 
measure. 



The maximum acceleration of 90% of all currently 
known pulsars falls in our acceleration search range, and 
as the 10% with higher maximum accelerations are still de- 
tectable during large parts of their orbit, where the compo- 
nent of their motion in our line of sight is less, we estimate 
P acc =0.95. 

As more than half of millisecond pulsars with known 
characteristic ages ap pear to be older than 5 Gyr 
l|Manchester et al]|2005l ) we assume recycled radio pulsars 
remain active over the cooling lifetime of white dwarfs. 



4 DISCUSSION 

There are several evolutionary scenarios thought to produce 
low-mass helium-core white dwarfs. In all these the progen- 



We will use our non-detections and the luminosity, 
acceleration and beaming fractions mentioned above to 
estimate the probability Pmsp that a low-mass helium- 
core white dwarf is made through the channel involving 
a millisecond-pulsar companion. We assume that our N=8 
candidates were correctly identified as low-mass helium 
white dwarfs. Taking radio frequency interference into ac- 
count we estimate that our pulsar search-algorithm success 
rate P eff is 80%. 
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In this case the assumption that our non-detection was 
the most probable outcome of our experiment: 

(1 - Plum X Pace X Pbeam X P e f f X P M Sp) N > j7 W 

leads us to conclude that Pmsp < 0.18 ± 0.05. 

Our upper limit for Pmsp is the first measure- 
ment of the probabilities in the formation of low-mass 
helium-core white dwarfs. In comparison, for general white 
dwarfs the mass-ratio distribution appears to peak around 
unity, strongly favorin g WD -WD over WD-NS binaries: 
iPinsonneault fe Stanekl (|2006l ) find that 50% of detached 
eclipsing binaries i n the Small Magella nic Cloud have very 
similar masses and lMaxted et al. (2002) show that WD-WD 
binary star s relatively often have m ass ratios close to 1. Sim- 
ulations bv lSchneider et all (|200ll ) predict the birth rate for 
general WD-WD binaries to be 100 times higher than for 
WD-NS binaries - a significant fraction of potential WD- 
NS binaries are disrupted in the supernova that creates the 
neutron star. 



5 CONCLUSIONS 

We have searched for radio pulsars around 8 low-mass 
white dwarf candidates from the SDSS and have found 
none down to a limit at the 5th percentile of the currently 
known millisecond-pulsar luminosity distribution. Low-mass 
helium-core white dwarfs can form around either a heavier 
white dwarf or a neutron star: we estimate the probabil- 
ity of formation through the latter evolution channel to be 
< 0.18 ±0.05. 
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